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ABSTRACT 


In a series of articles by Ditkin and Prudnikov, Meller, 

Osipov and others, an operational calculus for the operator DxD, and 
later for © xD: was developed and several interesting consequences 
and applications were pointed out. It also became apparent that the 
Bessel operator B = DxD plays the role of a derivative on the con- 
volution product 

x 

(fg)() =< f fGet)g(tar 
oO 


in the, sense that 
BCrse)"= (bijixe + t*(Be) -. 


This thesis starts by studying operational calculi for which 
BL = AxV, or in general B, = A(xty)V, play the role of differenti- 


ation. In each case a convolution * was constructed such that 
> * = * = . ° pee 
B, (fxg) = (B.f)*g + fx(B.g), i=1,2, . 


A transform L was also introduced in each case that satisfies 
b(ixe) = iiLe ~ and LB. f AUG, t0r “13 4.2. The firet operational 
calculus 18 essentially due to L. Berg but we aroved more results in 
that direction and included some novel applications. In case of B 


He 
and B, it was shown that any linear difference equation with poly- 
nomial coefficients is transformed under L to a linear differential 
equation with polynomial coefficients. This analogy provides a method 


of solving one of them by knowing the solution to the other. Moreover, 


in studying the first couple of operational calculi, we were led to 
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co n 
study transforms of the type Lf = ) = re] f(m) and 
Oo  (xtl)" 
n 
Lf = ) oars f(n), which are special cases of Jakimovski's 
Ca 


sequence to function transform, the [J,f(x)] transform. Operational 
calculi were developed. We also included some q finite difference 
analogues. 

Later in the thesis we studied expansions of powers of the 
discrete Bessel operator AxV; or more generally powers of any 

he x x 
eperator | tniat satisfies TC) =c(7.), ce =.0; ‘ec #0 for 
n n n- oO n 

n= 1,...:; in powers of the shift operator E and its inverse. 

We also introduce a discrete analogue of orthogonality with 


respect to convolution and used it to show that the zeros of the 


Rice polynomials are real and simple. 
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CHAPTER I 


INTRODUCTION 


Mikusinski developed his operational calculus via an algebraic 


approach. He embedded the integral domain (C[0,~),+,*), where 


x 
ti. 2) fxg(x) = f f(u)g(x-u)du_, 
Oo 


into its field of quotients. This quotient field is rich enough to 


x 
contain the integral operator f ECC) Cree | tee sG LO.) and hence 
oO 


its inverse, the differentiation operator In this setting, a 


‘dt . 
differential equation is nothing but an algebraic equation in the field 
of quotients. Solving this equation means finding a convolution 
quotient that satisfies it and if this convolution quotient is a func- 
tion, then it will be the solution of the differential equation. 
Mikusinski's approach has become a standard approach and 


several operational calculi were developed by using different con- 


volution products. Rajewski replaced (1.1) by 


x 
f f(u) g(x-u) du. 
fe) 


Qu 





fie?) f*g(x) = 


Qu 


x 


The resulting ring (C[0,~),+,*) contains a unit element, the constant 
function 1. It is clear that Mikusinski's ring does not have a unit 
element. Note that Rajewski's convolution also preserves dimensions. 


If we use the Carson-Laplace transform 


Glas) L[f£;x] = ae f ev t/x E(t) de... 
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then 
(1.4) L{f*g3x] = L[f;x]LEg;3x]_, 


wiere Lhe convolution in, (1.4) jis. that of (1.2). , Let ussdenote the 


d d 
operator an * ae by B. The operator B was called Bessel 


‘operator by Ditkin and Prudnikov. B acts on the convolution (1.2) 


as a differentiation. Indeed 
(1.3) B(f*xg) = £*(Bg)+(Bf)*g 


The obvious relation 


n n-1L 
a 


Sn es 


| * 


n 
= x : : Die 
shows that ar om the B calculus is the analogue of x in the 


1 
differential calculus. Ditkin and Prudnikov [13] introduced an 
operational claculus for the Bessel operator. This calculus enabled 


them to treat any linear equation in B with constant coefficients 


as an algebraic equation in a certain field. Their convolution is 


d une 1 
(1.6) GG) oe = 1) BCEN)GIATH Ey C1Mn) Teed 
ORF 0 


Under this convolution, we have 


xk ; ed . kts 


Pe a aa 








Meller [21] developed an operational calculus for the gen- 


=i) wit d 


x which contains 
dx rs ps aaa 


eralized Bessel operator BS ta 


Ditkin and Prudnikov's as the special case a= 0. 


i) 43)5. ° . < x 7 = (apet 43.4) 
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Amerbaev and Naurzbaev [2] introduced the notion of 
orthogonality with respect to convolution. The idea is to represent 
a known polynomial set POs orthogonal on [0,”°), as the Carson- 


Laplace transform of another polynomial set Q 6)» that is 


C197) Pd = ae if eae! 


oO 


x 
Q fede, x > 0 


Substituting for P fs) > from (1.7) in the orthogonality relation 


for Ps) we end up with a kind of orthogonality involving the Le Ss. 


Amerbaev and Naurzbaev [2] studied properties of the Q. 8 ofaC(lL 7) 


when Po are the Jacobi polynomials with argument (1-x)/2. Let 


608) 


e (x) be the corresponding polynomials. The convolution 


orthogonality relation is 





t-T ; 8B 
Oo tt Te 
(ees cf f ot SS ae Soar tay eee ee 
NS Ceedet pare [(8+1) Jie AL 4 2 it “2 


(x,)arjary} = patty (nese) g 


CG5r J eva, 
450 (t)Q, ni.(2ntatg)F(atBtnt+1) “m,n 


n 


They utilized (1.5) to deduce several properties of the polynomials 


(a, 8) (a,8) 
Q n 


= (x). Among other things they obtained 


{x)y irom those.of .P 


differential equations, recurrence relations and generating functions 


g6t>8) 


n 


for (x). These polynomials are not orthogonal, in the usual 
sense, on any subset of the real line and it is this fact that 
makes this new concept interesting. 


Berge [8] introduced an operational calculus on sequences 


using the convolution 


x 


“ograt 





n 


fxg(n) = ) () £(k)g(n-k) 
k=o 


and used the exponential generating function as his transform. Moore 
[26], Brand [10] and Traub [36] studied an operational calculus on 


sequences using the product 


£(k)g(n-k)  , 


Fh 
ws 
09 
COS 
3 
Ld 
Hi] 
ere~7 3 


which is a discrete analogue of Mikusinki's convolution (1.1). 
In [6] Berg developed an operational calculus for sequences 
based on the convolution 


n-l 


Lge } Sh te seo 


C128) fxg = 


oe l=) 


Berg's convolution in the discrete analogue of (1.2). 

In this thesis we give several discrete analogues of the 
Bessel operator and its generalization and construct operational 
calculi for them. 

In Chapter II we'consider the Berg operational calculus in 
more detail than was done before obtaining new results as well as 
indicating some novel applications. We shall find out here tHat 
@ discrete analocue for D x Dis An V. We then generalize Berg's 
calculus to a discrete operational calculus which corresponds to the 
generalized Bessel operator x” D ote D.j We note here that the 
operations which map the Bessel operators a the derivatives 


(analogue of Carson-Laplace transforms) are special cases of 


Jakimowski's sequence to function transform known as the; [J,£(x)] ~ 
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means [19]. In Chapter III we generalize the results of Chapter II 
by showing that a large class of Jakimowski type means induce cor- 
responding operational calculi by defining for each member of this 
class a convolution product, a Bessel-type operator and a Carson- 
Laplace-type transform preserving the basic properties of Chapter II. 

In Chapter IV we consider a discrete analogue of orthogo~ 
nality with respect to convolution. We uSe our concept to obtain the 
apparently unknown result that the zeros of the Rice polynomials [32] 
are real and simple. 

Finally in Chapter V we derive expansion formulas for linear 


operators IT satisfying 


x x 
Ge sn ee) 
( Disk @ © 
where Cos Cpe ces is an arbitrary sequence with he O. The 
discrete Bessel operator 8 is such a I operator with c., = j. 


Later, in Chapter V we develop, very briefly, some q-analogues of the 
(L,o) transforms and operational calculi. We also outline a possible 


way of extending the (L,a) transforms to several variables. 
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CHAPTER IL 


BERG'S OPERATIONAL CALCULUS 


2.1 It is well known that the action of the Bessel operator 


d d 
B=—t — jewski' i 
aE qe 08 Rajewski's convolution product 


ic 
feg(t) =< f f(u)g(t-u)du , 
O 


F F ee d A 
resembles the action of the differentiation operator a on ordinary 
products. The image of a convolution product under the Carson-Laplace 


transform 


-t/ 


L[f;x] = ec f air le F(t) dt ies Ones 


oO 


is the product of the transforms, that is 
kitegsx| = LPLix ile :x | 


d : 
Moreover L transforms B_ to ope a in the sense that 


| 
| 


Bie | = : Elet scc [oss 


and maps t Jn! to x", Thus the role of < and x” in the range 
of L is played by B and t'/n! in the domain of L. This 
observation led to some interesting results ([2], [13]), and so it 
is of interest to study a discrete analogue of this calculus. 

In the calculus of finite differences, the sequence 


{(.), j = 0,1,...} may be considered as the discrete analogue of 


x 
the functions { ae que QS. athe ites So Since 
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Cenk) Af) = fot). — £(n) 


Because of the relation 


ie: Le lees 
fie C4 ae 


we are interested in finding a finite difference operator 8, say, 


which takes os) to Mera) This is accomplished by 


e222) BY = A. EV. “ 
where 
CEs) te £CE) N= 6 (tb) —" £ (t=1) 


< k+l 
To see this we put, formally, B = ) a (t)A and use the require- 
O 


te 


ment 8.) ere) teres Oe). sc = COCSeE 


r~-1L 
k 
a(t) = 1] and a(t) =a =L eth, hark ee le lee tes 
so that 
ie k .k+2 he 
B=At+t d Gillette Weedgt Oo c6 


The operator 1+A is the shift operator E defined by 
(2.4) E. ECs) = f Cel) 


Clearly, 
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hence 


B. =A+tAV = Atv , 


which is (2.2). The assertion 

(2.5) B.C) Ses) Pe O;liece with ~<( 
can be easily verified. It is also easy to see that 
(276) B E(t) = “CerlyeCer l= (241) £(t) # tf(t-1) 


In the following we shall drop the subscript t from Ais Eas ve 
or Be whenever no confusion arises from doing so. 

We introduce a discrete analogue of the Carson-Laplace 
transform later. 

Let S be the set of all sequences Lene, with 
f(-1) = £(-2) =... = 0. S can be made into a complex vector space, 
in a natural way, by defining the addition of two elements of S by 
componentwise addition and scalar multiplication by af ={af(n) } 
for any complex number a. 

Throughout this thesis, for f ¢«S, by f(t) we mean the 
ce eomponent of, £ and formulas Tike (2.1), (.2)) and (2.3) should 
be interpreted as the sequence, in S, whose ae component is what 
appears on the right hand side. 

Now that we have found 8, the next step is to look for 


an appropriate convolution product that is well defined on S_ and 


satisfies 
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ose poe al tag, bad are eam 


ores B(ixg) = f*Bg + Bixg 


The convolution 


n 
(2.8) fxg(n) = V) £(k)g(n-k) , 


oO 
or equivalently 


n n-1L 
(2:9) fxg(n) = ) £(k)g(n-k) - } £(k)g(n-k-l) , 
O oO 


is a natural difference analogue to Rajewski's convolution 

d [E 

ae J £(u)g(t-u), hence is a candidate for the convolution sought. 
fa) 


Indeed we have 


Theorem 2.1. Ihe convolution product (2.8) satisfies (2.7). 


Proof. Clearly we have 


n-1 
fxg(n) = A ) £(k)g(n-k-1) 
e) 
Thus 
n-1 
(2.10)  (f*gtix8g)(n) = A} {g(n-k-1)B£(k)+£(k) Bg(n-k-1)} . 


[e) 


Substituting for B8f and 8g from (2.6) in the right hand side of 


(2.10) and simplifying the resulting expression we get 


n n-1 

(Gfxgt£*Bg)(n) = An{ ) f(k)g(n-k) - 2 ) £(k)g(n-k-1) 
Oo Oo 
n-2 


+ J} £(k)g(n-k-2)} 
Oo 


B(£*g) (n) 


oO 
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L. Berg [6], [7] has constructed, a la Mikusinski, an 
operational calculus based on the convolution (2.8). In the next 
section we mention few properties of that operational calculus and 
derive several properties that are not contained in [6], [7]. We 
shall refer to this particular operational calculus as Berg's 
Operational Calculus. For similar type of results see [14]. Note 
that discrete operational caculi based on different convolutions 


have, been studied in [8!},.[10}, [25],,[26] and,[27]. 


poeeere seOperational-Calculus. It is obvious that (8, +, *) .is 
a commutative ring without zero divisors. For brevity, we shall 
degotesthis vine by S also.esbet, {Y, t,.*}, or simply Y, be 


its field of quotients. The following formulas may be found in [7]. 


(ay td) PivOe0 sae ei fcn) beael Vein) 
cok?) {n}*{Af(n)} = {f(m)-£(0) } 
(2713) titi ieif a) }y=e1ofin) j., 
where 
Nn 
(2.14) ge(ny =e) E(k) 
Oo 
Let 
ee a 
(#245) Rone Gi; eee Oe levees 


The fundamental relation 


(2216) o Leet? = Sg E : Reyee OF lac 5 





1l 


also appears in [7]. Relations (2.11) and (2.13) provide representations 
for V and o in the field Y. The relation (2.16) expresses the 

fact that the functions aft) » a(t), a(t), ee» can be treated like 
powers 1, x, —— ~.+ » Formula (2.12) shows that A may be identified 
with the convolution quotient se on the manifold {f:f € S_ and 

£(0) = Of. Let e@{f(n)} = (nf(n)}.°- aft is easy to see that the 

operators 8 and 9 are not members of Y in the sense that there 


’ a ae a 
is no convolution quotient =-e« Y that satisifes o@f = qe mh 10Lr 


b 
pf => * £. 
Theorem 2.2. For f¢ 8S we have 
(2513) {n Vf£(n)} = ot, * BE - 
and 
(218) {nf(n-1)} = 01 * BE + a, xf 3 


Proof. By (2.13) we have 


n-1L n-1 
aereie a Y B£(k)} = { } AkvE(K)} = {nVE(n)} , 
oO @) 


proving (2.17). Now 


dy *BF = a, *a,*BE = a, *{n VE(n)} , 


1 


By (2317); so that 


n-1 n-1 n-2 
(a,*Bf)(n) = J kf£(k)-f(k-1)} =} kE(Ke)- )  (k+1) £(k) 
Oo O Oo 
n-2 
= (n-1)£(n-1) - } £(k) = nf(n-1) - (0,*£)(n) 
oO 


Therefore 


- : * 
he OTT = 
j el) + . ; 
4 \ : 
rs 
> y 
- ~ : 
= ' / 
- 
' p 
p 
4 hs 
- 
' ‘ fr 
4 i 
j “ 





“Feed. 3983- 





TZ 


(a, *B£ta, *£) (n) =n f{(n-l) , 
and (2.18) follows. 
(2.19) {nf(n)} = (a. +0.) *B£ + 0.4 ¥E , 


For the definition of convergence of convolution quotients 


we refer the reader to [7] and [16]. In this regard we have 


Sencar teers = (Ole eC) ero Ef. {1 51,1) cae ky 
n) 








then in order that lin exists, where ‘ SRA CA etry es, Le 
pda ghd Neato Saha yates Nandan, senor, aed 
n> © cee 
n-times 


is necessary and sufficient that | £(0) | ol Moreover, it this 





dimit exists, it willbe zero. 
The proof will be postponed until the end of the next 


section because it uses some results to be developed there. 


2.3 A Discrete Analogue of Carson-Laplace Transform. Now we come to 
the problem of finding a discrete analogue of the Carson-Laplace 
transform. Denote the required transform by L. We require L to 

be linear, to map S into the ring of formal power series and to have 
a nontrivial range. We write L[f;x] to mean the image under L 


of f(«S). L must also satisfy 


(2.20) 7 Dee ‘ eres) 
and 
(2.21) Lifeess) = Litsxl Lexis 


The linearity of L implies that it must be of the form 






; im 
. week Lee : = 
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7 “5 is an = a 
1¥ f Ad? } 1 

- i 
1 ‘ 

ys 
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(20225 L[f;x] = } a Cait), 
Oo 
where a(x), a, (x), ee. are power series in x, 


Theorem 2.4. The transform L is given by 


© n 
x 
(2.23) L{f;x] = } erap ile 
o (1+x) 
up to a shift in x. 
Proof. “Set 
6, Cx) = Lf, 3x] ; ee =O LL Soro nats 


The relation Oo Pasere erorestb afecss. implies. by (2.21), 
Lo sxIL[f3x] = 1s fee) ees 
and hence 
85 = Lfa 3x] =1. 


Using (2.20) we get 6, (x) = 6 (2) > that is 


6, (x) = xtC, C -is constant. 
Recall that a = O14 XO Feo KO) PcOGieic 0 eee Thuct by) (2.2 1) and 
Ne ——$~_-—_—_——" 
k-times 


the definition of 6, Cx) we have 


6, (x) = {Lla, 3x1} = (xtc) § x 


Thus 


‘ 
: f 
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r= 
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ole eae 
(xe) = & a) G) 


Multiply both sides of the above equation by (-1) 57" (5) and 
n 


summing over all j we get 


n foo) 9° “ } 
Perse (2chC) Meera a(x) J «qn? 


(2.24) 
Ciena k=o j=o 


and by virtue of the well-known identity 


Sey oe Aone 
bo SENS A ecb 


this formula reduces to 


a Gx) = (xtc)? /(1txte) ttt ; 


and we take c= 0 since, in general, we allow a shift in x. Now 
it only remains to show that (2.20) and (2.21) are satisfied for 


Pe Bee (S's 

a n n-1 , 
SS See £(k) g(n-k) - €(K) g(a-k-D)} 
n=o Clits) sek Ue ; 


L[f*g;x] 


a x ak.» x 
} 10h Gua, L eee 


n 
x g(n) 
5 rap Sea 


a eee le 
a ve ORG ) 


L[£;x]L[g;x] 


Finally we have 


boy (~) (i~) fveolg 
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(1+x)~ 


- f(n) {n + ae 2n-1 + (nt1l)x/(x+1) } 


pen 
Ont A.-8 


& 


ii 


(a)? J GE" e@) @- BS) - Stes , 


and the proof is complete. 

From now on, we adopt (2.23) as the definition of the trans- 
form L. We shall always use small letters to denote members of § 
and capital letters to denote their transforms under L, that is 
£G:penwill stand for LEGS). 

The above L transform appears in Berg [7] and Ditkin and 
Prudnikov [14] with a change of variable. However, no attempt was 
made there to derive (2.23) as we did in Theorem 2.4, 

Now we look at transforms of the fundamental finite 


difference operators. The formulas 


(2.25) Tint(ny x] ie GEC DB Whete , D‘= = 
(2.26) ; L{VE(n) 3x] = F(x) /Cetl) , 
(2.27) LPAP sx] = (F(x) — FO)}/x | 


and 
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C2220) Rlotsx| = Gtr) 5 


follow easily from the definition of L, hence we omit their proofs. 


By induction we get 


(2.29) PIHRECn Youd oe Tentetl Se GO oF 
(2.30) Byer Cay te hes PO LCS 
and 


k= = 
kr Lg(-D) 


Gees) wih cet has CEOs) 8 (0) Se CO)-Y oes Une > On. 


wnere (OQ), £'(0), ..., are defined by 


@) 
The formulas 
(2.32) a es 
and 
(2.33) Cee ieee £ Gt Die Fr) 


follow also easily by induction. 


The formula 


nf ok k 


(2534) it £3) 7x p* Pex) eee 


where 


(ie) . T for’ k®"0 
> 
n(n-1)...(n-k+1) for = Kkeo.0 


may be proved as follows. k= 1 follows by (2.32). If (2.34) is 
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valid for some k, then 


kt+1) k+1 Ke) ok (kh) k 
bin’ dy Pseloe L[ (n-k)n yg i ae xD(1+x) Ln y*ve;x] 


- xin’ pk vex) 


NCU Df Rt) Cae) foe DEF Ca) (lds) 


eee Tn ery Hees Cl SFG) [ck 


k+]_k+ 
_ x 1k 8 


(x1) (F(x) /Cxt1)} , 


and the induction is complete. 

Relaciones (2.29), (2.40) and (2.31) show that the L 
transform of a linear difference equation with polynomial coefficients 
is a linear differential equation with polynomial coefficients too. 

In particular, a linear equation with constant coefficients is 
transformed under L to an algebraic equation. Relations (2.32) 


and (2.33) show that difference equations of the types 


|] P12 


g 
a, { (at) V} f(n) = g(n) and ) a Cay ea) = e(n) , 
k=o k=o 


are transformed to 


Ke 
a { (tL) DEF (x) = G(x) and ) a, (xD) “F(x) =WiG( Xi 
k=o k=o 


il D1 


getty) are constants. Furthermore 


respectively, when aye a Q 


ahs 

according to (2.34) the difference equation 
Q 

k) _k 

y an yey = g(a), 


k=o 


is transformed to Euler's equation 
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Q 
) ay _ K(x) = G(x) 
=o 


The above correspondence provides a possibly useful method 
of solving difference equations because the resulting differential 
equation is very often easier to solve than the original difference 
equation. 


tet X*= {f:f_< S. and sup |£(n) |2/" < 


o} . Thus, for 
£f e X, F(x) will be holomorphic in a neighborhood of the origin. 


t is plain that X is closed under addition. X- is also closed 


under convolution. For, if n> 0 then 


n 
| ) E(k) e(a-ky ee Sten: sup digesee kee Oe) the 
k=o 
x sup Abas SSP eM As 
ee 
where 
Pek 1/} 
A = max {1, sup |f(k)| / , sup |g(k)| ie £(0) (0). 
k>o k>o 
1/n 
Therefore sup | £*g(n) | < © and fxg ¢ X. 


n 


For f € X, or equivalently F(x) is holomorphic in a 


neighborhood of the origin, the inversion formula 


(2555) f(n) = a f eee ere Gane: dz 5 
c 


or 
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es Se eA yyy 
(2.36) COME cm ye atl (l+z) dz , 


with an appropriate contour c, follow by Cauchy's theorem. Note 
Pig tgte met ce etien sOnwald. be, Aft ect. Ve, int(n)}, =... and all 
the operations we performed so far, for example rearranging terms 

in series, are justifiable from the theory of functions point of view. 
That L is one-to-one follows from the identity theorem for power 


series. Note that the proof of 


(2.37) LI) 3x] = , 


is contained in Theorem 2.4. Relation (2.37) provides an easy 
co 
k ‘ 

inversion formula. If F{x) = ) ax , is any formal power series, 

n * fe) 
then f(n) = ) aC). 

Oo 

It is also easy to see that 


foe] 


(2.38) Ulex) = > x A £(0) 


ie) 


In solving difference equations one might use the operational 
calculus approach or the L transform technique. We would rather 
use the L-transform in our illustrations and examples. 


Now we are in a position to prove Theorem 2.3. 


Proof of Theorem 2.3. Clearly 
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<0) 


Therefore (0) = {£(0) }* for all. ks ‘and ¢) (gy converges, 
as k+>™, if and only if |£(0)| <1 or £(0) = 1. In the latter 


case, £(0) = 1, we get 
6) (1) = THEAE BOYS EQOVIS = 1+kA£(0) . 


This implies Af(0) = 0 and by induction we may deduce that 
AX(0) =) for *k > 0. Therefore; in this case, £ = {1} which is 
a contradiction. 

| Conversely, let ee(O) sale then (2,39) , for “k/son, 


implies 


ais ) 
eth tt 


j j 
Hore < |£(0) | [apt eco) esa DEO), 


(k) 


which shows that f (n) converges to zero, as k++, and the 


sufficency part follows. 


2.4 Examples and Applications. In our first example we solve the 
three term recurrence relation satisfied by the Hermite'’s polynomials 


HG) 
Bg ao eve) Zo FC). pO elas 


with H_,(y) = 0, H (y) =1. The substitution Hj(y) = 6 (y) (2y)” 
leads to 


2y" ho. = “BG ky : 


Applying the L-transform to the above equation we get the differential 


equation 
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O'(x) + O(x) (2y“tx) J 2 2y7 7x2", 
where 
O(x) = L[0 3x] 


The solution of the above differential equation is 


poe x Ze 
xO(x)e ” os Dy” f er arenes 
fe) u 


2 


; 2 
If we perform the change of variables vie = ay + t we will 
obtain 


x0(x) = y f ara: + vfs 
O 


Se?) ae 


and by repeated integration by parts, we obtain 
a 1 P Peet eA 
GG) eats Ga ls tst-x (ye 


so that 
aie tah 2 ye 
ay) 2) 5G) ene er 
j=o 
and 


=9 
emt Coe TY 


N[5 
‘ 
ws 
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Ne 


H G)== (2y)" GF 


which is a standard hypergeometric form of HG). 
Our second example is to solve the three term recurrence 


relation 


(2.40) (nt1)y_,, Cv) = (2ntl-uta)y (u)-(nta)y _, Cu), bn l, Ae, War 0, 
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whose solution is known to be the Laguere polynomial Vaan The 


recurrence relation (2.40) may be written as 
By Cy) = ~uy tovy “ 
whose L transform is 
Y'(x) = —uY(x) + a¥(x)/(xtl1) . 
Therefore 
Vex seClixl. 6) aew 


which is the standard generating function 


; PS a2 r a-1) 71 a xt/i-t 


As another example we solve the three term recurrence 
relation 
(2.41) (ntL)y_, (uw) = (2ntL)uy(u)-ny,_j(u), 
Cao eevee end wie, Cu)ieaa bas 


The solution is known to be the ae Legendre's polynomial 


P ow: Rewrite (2.41) in the form 
= = 
By = (2ntl)(u-Dy, 


and apply the L~transform to get 


(1+2x) (u-1) 


Y'(x) = 5) 
1+(1-u) (2x+2x) 


Y(50) mae 


Therefore 





L/2 
(2.42) Y¥(x) = {1+(1-u) (2xt+2x Ags BA ) : 
which is nothing but the generating function 


) P (x)t" = (1- peti: ) ~(1/2) 


Furthermore, by expanding the right hand side of (2.42) in powers of 


x we obtain 
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One can derive another hypergeometric form for the Legendre 
polynomials by using the substitution y 6) = uf (u). The f's 
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Thus 
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2 2 : 
F(x) = {I-x(u-D/u"} 2 = J aE eX , 
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Our next application is to solve the non homogeneous linear 
difference equation with constant coefficients using the L_ transform. 
The uSual method of solving this type of difference equation relies on 
solving the characteristic equation as well as obtaining a particular 
solution. However, solving an algebraic equation or expressing its 
roots in terms of the coefficints is indeed very difficult in general. 
Our approach avoids this difficulty and we obtain the general solution 
expressed in terms of the coefficients in the difference equation. 
Let 
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Let 
k : il co : 
(2.45) {1+) (-l)°a,x } = J ax 
L : fe) J 
From (2.44) we get 
k=1 k-ie1 = 
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F(x) = {x G(x) + J Ci) a, ) ma Fe) (Oy) oe 
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It is clear that the above solution contains k arbitrary constants, 


F(0), F'(0),...,FOe 


(0), which is the right number of constants. 
The first term in the right hand side of (2.46) is a particular 


solution while the second term is general solution of the homogeneous 


equation. Note that the h's in (2.45) may be given explicitly as 
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follows from the multinomial theorem. 


As an example, consider the difference equation 
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(2.47) Cale =!CGo.+ Cachan, 


that was studied by Weinshenk and Hoggatt [37]. They solved it by 
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an operator method as well as polynomial expansions but both solutions 


are rather lengthy. We first write (2.47) in the form 


2 r 
NER pale DAY Shp stilted 10 one 
nee, D 


So that a. ”°="qa,"="=1.°" It is well now that 


n 2, =1 
FP = Clix x) : 
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where Fo oFio--- are the Fibonacci Numbers. Thus the Bs of 


(2.45) are given by s = Gone 
at j 
ark KB ae. 
j (1) ae 


and the solution (2.47) is 
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since FY and FA are linearly independent solutions of the 


cell 
corresponding homogeneous equation. 
Our final application of the L transform, in this 


chapter, is to Bateman Polynomials F(z). Bateman [4], [5] studied 


the polynomials 
= = r= ‘ » ° 
F(Z) 3 Fy Kekbhl (zt )/2c bee lye, 


quite extensively. In what follows we use the L_ transform to 
derive some of his results as consequences of well known properties 


of the Legendre polynomials P(x). Recall 
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Let n = -—(1+z)/2 in Bz) and denote the resulting function by 


f.(m). Clearly 
(2.48) L[f, (n) 5x] = oF, (-k, kt1315-x) = PY (2xt1) 


Relation (2.48) enables us to derive several properties of F(a). 
Under vig Legendre's differential equation 
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n 2 
J4n)*xB £ Me 7S = ewig P\ es 
((,)+n)*6 £ (n) + (it2n)*8 £, (n) k(k+1) £, (n) 
which leads to, after some manipulations, 


(2.49) 4 (2-1)°F (2-2) ss + (241) °F, (242) : & {(z-1)#(21) FB, (2) 


= k(k+1) Fy (z) 


So, (2.49) is valid for negative odd integeral values of z, hence 
FCueALIE 2. 
Using this method we can also derive the generating 


function 


“wan Ball ae Pa 
EF (G> 5) + 5 are 


| oe —1 
Pea) ete 


from Bateman's tF (z)}° for the standard generating function 
) Bac) oe = (i-2xtee?) 7! A 


for the Legendre polynomials. 
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iN 
The general Lucas polynomials tee (ajs+++sa), or 
simply oe are defined by 
CN) Reale ate (N) nm. 2() 
NERS OF hese c ooeaptey og Mi gaa pemeaa te ; 
Gayo ee Cte | glee (8) Nek 7 
with U5 = uy SP hover UN? Ow, Uy-1 t.vy Lhey arose,in 


connection with finding the wae power of an.) Ny xeNoematrix: [3]. 
Barakat and Baumann [3] expressed a need for a closed form which 
was obtained in [30]. The required closed form obviously follows 


from (2.46). 
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CHAPTER IIL 


EXTENS [ONS 


We note that the finite difference analogue of the Carson- 
Laplace transform, that was studied in Chapter II, is a sequence to 
function transform known in summability as Abel transform. Since 
this is only a special case of a much more general sequence to 
function transform, the so called [J,¢(x)] transform introduced 
in [19], then it is of interest to see which, if any, of these 
general transforms yield an operational calculus, or a transform 
theory, along the same lines we developed in Chapter {I. We shall 
restrict ourselves to such $(x)'s that are power series in x. To 


be exact let 








_ + (-x)™ ce 
(3.1) Dea ie ee Oo et (a) 
) ‘ dx 
where 
(3.2) Se ean 
Oo 


Substituting for.-$6(x) in (3.1) and using 


a ea oe kane 
(3.3) CG me ee ete (0)... 
oO 
we get 
(3.4) L{f30,x] = J} (-x)" a, a" £(0) , 
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provided that we can rearrange the terms in the series appearing in 
(3.1). At any rate for a given sequence a: Crt sreesG seer we 
shall take (3.4) as the definition of our transform instead of (3.1). 
The dependence on the sequence fo}. is explicit in the notation 
L{f;0,x]. If (x) has a positive radius of convergence then (3.1) 
and (3.4) will be equivalent. 

In the case $(x) = Gees or rather nS (ely we 
have the L transform of Chapter II. Whenever we drop the dependence 
on ae that is we write L[f;x], we shall always mean 
LL£3;(-1)",x] of Chapter II. 


We note that formally (3.1) is equivalent to 
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Lifso,x) = J (-x)"t@ J Oya. , 
n=o j=0 


which in turn has the operational representation 
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from which we see that 
L{f;0,x] = F(~-xE) Ot 


where FC )e= Cf *x). 

As before, we shall use small letters. f, g, ... to denote 
sequences and use capital letters F, G, ... to denote their transforms. 
We shall refer to these transforms as (L,a) transforms in general 
or (L, a.) if we are spedifying the sequence (one : 


It is natural to require (L,a) to satisfy 
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which is the case if and only if oe 1. We also require (L,a) 

to be one-to-one. Relation (3.3) shows that f£ is uniquely 

determined by £(0), AEC), ..., hence our transform is one-to- 
1 


one if and only if none of the a, 8 vanish. That is why we 
I 


shall assume 


(325) a =1 and Cer teatime 2 ee te 
fe) n 


Now we look for a suitable convolution product and a finite 


difference operator BO such. that 


- - en ; 
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Recall that we referred, and shall always refer, to the case where 


= (ists as 8. 


3.1 Properties of the (L,a) Transform. We define the sum of two 


sequences as before. The requirement 
Coro) Piers. <i neer Cx) GX). s 
defines the convolution product uniquely. In fact we have 


Theorem 3.1 The relation (3.8) is satisfied for all 
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e n 5 et 1 roi 
(3.9) Ere) ene) 2 Meer (O) A, e 3 (0) 
Proof. [f (3°38) is satisfied, ‘then 


(-x)"o a" (e%g)(0) = { J a, (-)*a*e(0)H J a,(-x)Jadgcoy} 
i=o j=o J 
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so that 
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; al 
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aL n-i 
; a, iA ECO)A AGED EO 


auaeto. 9) ToOllews by (3.3) 

Conversely, if (3.9) is satisfied, then (3.8) will follow 
by elementary series manipulations. 

So, we adopt (3.9) as the definition of the convolution 


product. Let us define a sequence CooCyores by 
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and introduce two linear operators AY and BA by 


n+l 


ag ia oat 
(3.12) Seas etn) se on Ce tiee ome (Or, 


and 


(3243) BY = AYR Va, thie, (8 £) (n) = A nv£(n) , 


; ’ wr r 
n T ’ bi 4 + . 
+: e One =f : a 
+ A. 
h 
- 5 a ee ee 
5 F 
zs i t | 
2 : : a re = 
I ) = 
, _ i i = 
. tOxe wt 4 ey ae 
: y 
2° i , a ¥ 
‘ an 0, 
» ; 
res ff ‘ | fic ht t) 11 
1 
A) 4 
iA « mgt } i 
rth} . ) i \ h fi ‘ ‘4 J te) Ue 
ae fi ; : 
: (i : . OF ~ 
Pi = c -* 
—* % s 
= : = : 
fr 
ieF 
64 4 
a 
i 
at 3 ee “ay s* i a, J , i * 
» a - © 
: oe] Go 
i. 2 {* “are ‘ 
= ~ L 
r > 2 u! ts ae 
0 
; . = i Ta tad 
' we « 4 en . | 5 eins s = 
g - * . sg 
y 
= Ps 
‘ - Me DPE v 
i> i . 
* 
. . b 





‘ ‘ : - 
@ i - him. * A> ejoteeagqe genntl ori sockorss 
‘ 


The operator BY satisfies 
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(3.14) Geta = (REI) c 
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n KtL 
k+l. C 4 £(0) 


To see this, we have 


n 
Byf(a) = Alnve(n)] =) Gop ho [m(E(m)-£(m-L)) 1s 


and by the finite difference analogue of Leibnitz formula, see 


C341, p-: 11), we obtain 
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= 2 Goa A FO) 


proving (3.14). 
iteis ciear, from (3:12) and (3.14), that for the sequence 
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It will turn out the By as defined by (3.13) is the sought operator 
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(3.20) L[A f ho oc] = FOx)-F(0)_ : 
Q xe 


as well as (3.6) hold. 


Proot. (3.1/7) follows from ree TP Seg, anid s(o.4) 
——— rc 'n=0 r,m 


By (3.13) and (3.3) we get 


m : mt) Om mbl 
A A ECO) c+ £(0) = - Pie A £(0) 


Thus 
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L[A f30,x] = i (-1) Ot 
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peeviue (3.13). *Gimilarty (376) tollows from (3.14) and “(3.3): 


Corollary. BY Satisties (3.7). 


Proof. This follows by (3.6) and (3.8). It can be proved 
divectly iirom 4(3. 79) eands(3s14). 

Note that if c. ARSE | po lyneniaiein n, then AY will 
be a finite difference operator, for example ce i Dh Sry: 
correspond to A ae AeA and yA + AnV respectively. Jn fact 


we have 


Tnaecrem 3.5. = Ik i P(n), where P is a polynomial, 


then A, a oP (iA and BY = P(r) 6 where tT = 1+xV. 
The proof is rather easy and is omitted. These operational. 
formulas can be extended formally to functions which are power series. 
At this stage we can use the (L,a) transform to solve 


equations in AY of the type 
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with constant coefficients do» b> tsen0De. Let usssielus trate the 
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method and its advantages by an example. 


Example. Consider the equation 
(3.19) AnVf£(n) - f(n) = g(n) 
We take a = AnV , that is c, =m, so that 


F(x) = L[£;(-1)" nt,x] = n! x A" £(0) 


mi 8 


n=O 


Applying the Gately an) ) transform to (3.19) and using (3.18) we 


get 





£(0) # xG(x) 
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F(x) = 


Thus, using (3.17) we obtain 


nh n 
£GD = £10) Gb /kKi+t J. Coes(n)/k! 
k=o k=] 
The observation 
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eee) () /! = (1-8, )/4 ; 
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leads to 
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On the other hand if we use a generating function method, i.e., if 
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we put 8(t) = ) £(n)t", to solve (3.19), we will get 
Oo 
L/(1-t is 
Diner et! ( : j ere isu) any du_ 
(1=t) : I-u ° 
) 
where 7n(t) = ) g(n)t. Now it is apparent that recovering f(n) 
Oo 


from @(t) involves very messy expansions. The use of an exponential 
generating function leads to a second order differential equation. 


The use of the L transform of Chapter II leads to 


s 18 
Poe abbeseldt 
e 


s x 
L[f;x] =e 
and recovering f(n) will again involve messy expansions. 
Now we go back to the general (L,o) transforms and look 
at the transform of nf(n) and Af(m). First, we introduce the 


linear operators C and C, defined on formal power series by 


n cae 


(3.20) Cx =x ple ba Ae Bo Got eh 


(ntl)c 4) x : Tes Ole ome. 


(3.22) lege 


Theorem 3.4. Under general (L,a) transforms we have 


(3222) ULPAL+ a,x) = CeLLesasx)) _ , 
and 
(3-23) L[nf(n) 3;a,x] = «£1 [ f50,x] tex COLtt ex) a 


Proof. We have 
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' % . eh ag bp 8 reat atl n+]. 
L[Af;a,x] =) (-x) a A" £(0) Y (-1) ROL niece ag Mmr(0) 

) ) 

nth +1 + 

=) G17 a,a7 ecoyctx +] = ccuftsa,x]) , 
i nt+1 
proving (3.22). Now 
m ‘t m—-1 
A nf(n) = mA Cl youee 


n=o 


follows from the finite difference analogue of Leibnitz formula 


([34], p. 11). Therefore, since m A™+£(1).= ma™e(0) + m™*£(0) , 


L[n£(n) 30,x] = ) (-x)"a_{na™e(0)+na™”"£(0)} 
oO 


a Kou a hag 1 n 
ig ese { Y (-x) aA a Cae: 


° ntl n 
: , (Se) ee(5c) s(t) a £0) 


the right@hand side of (3.23). 


Teo Usup {Jaa"e¢0) [1/ ee=al oe Aor othe function 


F(x) will be holomorphic in some neighborhood of the origin. Let 


vs Reis : sup jo a™=¢0)|7/™ < 


m>o 


orn. 


Y is closed under both addition and the convolution product (3.9). 


For, by (3.9) we get for m> 0 


+, “) oe res 4 . oo a EA, aes iain me +e aap! ’ 
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a A™exe(0)|1/™ < Gly anak {|a,4"£(0) | 


x max {{a,07¢0) | 


where v = max {1,£(0),2¢(0), sup Jo,*eo)[*/*, 
i>o 
Thus f*g € Y for heeetveyandy see ry. 
By using the argument in Theorem 3.1, 
fore fic°Y} awherée, in thisssetting}sequalities 


etc. must be interpretted as equalities between 


1/m 


1/m 
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T10es jj acent 


[A 


sup Ja,atecoy|/93, 


j>o 


we can prove (3.8), 
as 63.46) > 338) r) 


holomorphic functions. 


With this in mind, Theorem 3.2 is valid if the left hand side of 


(3.6), (3.17) and (3.18) are holomorphic in a neighborhood of the 


he poe ps nee 2 mn 
Origin.  Inis is indeed so since A (eo) owe =6 


m 
aA A £() = pe er 


imply Toy % =n0ebeat -joeaYn, andathat iboth ALE and 


members of Y whenever f is so. 


> 


r,m 


mt1 m ~ mtl 
A £(0), and aA B FC) = ~Crtl)a, 44 £(0) 


8 £ are 
a 


This approach is undoubtedly more restrictive than our former 


one, using formal power series. 
allows us to use complex function theory. 


formula 


F(z) f 


; Et) dz 
271ia 
n 


n mq — 
(3.24) A’ £(0) = ove 
i 2 


or equivalently 


On the other hand the present approach 


For example the inversion 
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25) flay aa hy 
¢€ 


n 
: kn, -k-1 
one Lb (G2 fa.) Fda , 


k=o 


is obvicus, with F(z) holomorphic in and on the contour c. 
If we restrict ourselves to operators AN or Sequences 


co ea 
fob.) formwitche, Of<tsun ales sac me cnen both Ck(x) ) and 


no 

n7o 
CF(x) are well defined and holomorphic in a neighborhood of the 
Origin for f ¢ Y. Moreover, Y will be closed under finite 
differences ard products of its members by polynomials. In this 


case (3.22) and (3.23) follow by a legitimate rearrangement of 


terms in the Taylor series of the respective right hand sides. 


3.2. An Operational Calculus Approach. In this section we develop 
an operational calculus based on the convolution product (3.9). Our 
approach is, as before, an algebraic approach, see [16]. It is easy 
to see that the set S under addition and the convolution product 
(3.9) forms an integral domain. We shall denote its field of 


quotients by Q. 


Theorem 3.5. We have 


n 2 —_ 
(3.26) A E-f) =f, 
if 
and 
(3.27) (R-}xa £ ane HR RG) SOO 


a: 





Proof, It is easy to see that A £(0) = 6,4, if f(n) =n, 
Thus 
nN a a 
A,CSo}ef) = a ty Cy Et to! 6 (0)} = £ 
10) (&, r S a 
if r=1 Leas 


Relation (3.27) can be proved similarly. 
Relation (3.27) shows that A can be identified with the 


quotient {1}/{n} on sequences {£(n) }° with §£(0))= 0, 


Definition. The index of a sequence as denoted by i(s,) 
or i(s), the index of its first non-zero term, that is 
i eS i 3 ) 0 ° 
(s,) min in s, # O} 
‘ ; ; ; i 
it.is clear. that i(a) =a min {i : A a, 7 O} . Theuseful 


relation 
C3828) itatb)>: Sita) 1th) ~ ; 


follows from 


5 


r < 5 8 r-i 
aA £xg (0) ) aja, £(0)A “g(0) , 


ie) 


which in turn follows from (3.9). 


: d a 
Theorem 3.0 In order that a convolution guotient b 


belong to S, it is necessary and sufficient that ia) > ib). 


Proot WL: cr = é€ S, then b*c = a and so 


o|m 


i(b)+i(c) = i(a), and the necessity follows. 


Conversely, if i(a) > i(b), then the system of equations 
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a A - v ) a, A x JA D. > i Se OL eA sens 


f n 
in the unknowns Xo» AW AP trey iter 
Oo 


. » »«» has a unique solution, 


proving the sufficiency. 


(n) 


As before, we denote Peigea, eee Oy £ 5 
———— 


n--times 


(n) 


Theorem 3.7. Let {is te Seewith 5 Pods Wothenve linet 
no 


exists if and only if Is, | — 1, in wiich case the limit: will be 
Zero. 
The proof of this result is similar to the proof of 


Theorem 2.3 and is omitted. 


1.3 Special Cases. In this section we study two special (L,o) 
transforms. Recall that the transform L of Chapter II is another 


special case. For more applications see Chapter IV. 


1.3.1 The Case. a (-1)"/n!. In this case nes Liriee eGo >. 0, 
or o(x) = e*, This special transform, L[£;(-1)"/n!,x] may be 


defined by 


2 nN 
Dles@elya/ntxlem es 4) f(n)o., 
QO 


and is indeed a modified exponential generating function. This 


transform, apart from the factor e 7 was studied by Berge [8]. In 


this case the convolution product (3.9) reduces to 


fxg(n) = 
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()a*g(0)e(n-r) 
oO 


irs 


abt 
tA ® 
met 


fe 


ie eave. , 


pe eg FY 


ae 


a 


J 


A 


- f 
Bt Apaogee. | 


’ 


4 


way i hi Ort 


7 
rh 
ro 
© 
‘ 
we 
te* 


. ae 
swoniny $69 \ak 


i 
- 
. 
at Ge 
_ 
a 


m4 
oo 

; 
< y 
f 
bel 
So 


a hodE Yas 
" w 


a 









* - 
- 


é 








the operator C is the differentiation operator and C is the 


identity. Thus we have 








n n 
(3.29) De ee en ee] nn 
Ne The 
and 
n n n 
(3.30) Liné (n) Sx] : xi [63x] “ xp[ és? et 


For more properties and numerous applications the reader is 
referred to [8]. 


As an illustration let us solve the recurrence 


i303 1) E42 = (nt4)£ Gre, (re Oe with fo =a and fy =b 


which appeared in the advanced problem section, Monthly, vel 80 (1973), 
problem 5911. 


Fitst we rewrite (3.31) as 
A-£. = nAf +2Af£ +2£ +(b-4a)6 : TieccenO-gebvoreie« 
n n me Do n,o 


By (3.29) and (3.30) we obtain, for F(x) = Lik-l)pin!.£.x) the 


differential equation 

(3.32) (1-x) F" (x) - (24x) F! (x)-2F(x) = (b-4a)e™ 

Multiplying (3.32) by e* and integrating the result we get 
eX (1x) F! (x) -20*F (x) = (b-4a)x + constant , 


and the initial conditions imply 


(3.33) (1-x)F!' (x)-2F (x). = (b-4a)xe “tae ~ . 
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The solution of (3733) is 











5 
2 es 
Cl=x) FG) = { {(b-4a)ttale Aaiank + constant . 
Oo 
Therefore we get 
aoe Z (b-4a) 
F(x) = 7 L(b-4a)x"+-(b-3a)xt(b-4a)} - 
(1-x) (T=x) 
and 
° “il 
f= 4a-bin(a!){3b-llatGa-b) } Gt) 7} 
fo) 
As an application of this transform we derive Kummer's 
transformation 
= x a ® a 
pF, C3Psx) e pF a3b;-x) 
We write 
S Teh © (a, )(-n) n 
pPy Ca3bs =) = ) ey = Bl ) ey ee n] 
‘ x & oe Ke ine 
(2 ye -x ¥ ae n 
iD Strom £4) Re rere Ps 
(ba) | ni 2 n (b-a) | 


Y+n 


1.3.2. The Special Case. Ot = ibe \ Boel thie .case 


o+n - : a 
Se el corny vee 0. We also assume that a is not a negative integer. 


We shall denote the operators AY and oo in Pihis*¥ease-by ay and 


oe and shall not use the symbol y anywhere else. Clearly 
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f (n+1)-£ (0) 


(3.34) ae? = Af(m) + y a 
and 
A205) ee = BGs) 2). 


It is clear that our B is the finite difference analogue of the 


-yd_jytld 


alized Bessel : eet 
eeneralize essel operator at Ae 


introduced by Meller [21], 
since the latter can be expressed as Sy te “). Note that in 


case y = 0, the operator ) reduces to A and hence our operational 


calculus and transforms reduce to those of Chapter II. 


ie 


For simplicity we use Tel to denote LI[£;(- te = ( eat: 
Clearly 
i206) J Ifsx] = * xn C1") A"E (0) 
is r 
and 
+ 

; x" Care) 

Cs237) J7Pe sx] es ————— f(n) . 
nF fe) (ig) Gao 


Relation (3.37) shows that the a transform corresponds to the A, 
method of summation, a generalization of the Abel means, introduced 


in [9]. By straightforward manipulations, we obtain 


(3235) a [V£(n) 4 oa eH x] = TS la) ‘ 


Relation (3.38) is a characteristic property of the ve 


transform. Indeed we have 
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Theorem 3.8. If an (L,a) transform satisfies 


LIvé(x) + —& £(m-1)30_,x] = h(x)L[£30,x] 
nty n 


whenever L[f3a,x] and h(x) are holomorphic in a neighborhood of 





neat ‘ . : = 
the origin and yy is not a negative integer, then h(x) = (1tx) 


aude.) ese), sup ce a change of variable x=» )x. 


Proof Leis cl camthateueLif:o.x}) willobe holomorphic in 
a neighborhood of the origin if f(m) vanishes eventually. For 


those sequences we have, by (3.3) 
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identity we get 


ntyt+1 


ve aie aus Pate se sary ss 
Therefore 
ney eh 
Sagi t tty oat pane 
hence hy #0 and we may assume hy = -l1 since we allow the change 


of variable x-> \x. Consequently the (L,a) transform is nothing 


but the i transform. That h(x) = (lyr follows from (3.38). 


Corollary. The only (L,a) transform that satisfies 
Live coves hts) Lt vosx) ek 


where L[f3;a0,x] and h(x) are holomorphic in a neighborhood of the 
origin is the L transform of Chapter ITI. 
As an application of this special transform, we use it to 


derive the formula 


Zz 


= = mic ee SO 35) este 
£3539) oF, Casbse3z2) = (l-z) oF, Case Doce i 


Let £ be defined by 


(ity), (©), k 


C3.40) J Lisx} = igs (-x) 


oO 


Therefore 


2  (b), (=n), 


f(n) = ) k! (c), 


na = (cb) /(c) A 
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en nds 
CY), 


‘ coal § 
since Jy maps x to 


Chie: (c-b) x" 


o (on! (c) Geo ee 


| 73 8 


(3.41) J (£3) = P 


“Yoh elit Xe 
(1+x) oF, Cytts¢ bees ree 


Relation: (3.39) now follows. from (3.40) and. (3.41). 





CHAPTER IV 


ORTHOGONALITY WITH RESPECT TO CONVOLUTION 


In this chapter we introduce and investigate a discrete 
analogue of the concept of "Orthogonality with respect to convolution 
introduced in [2]. We also point out discrete analogues of Rodrigues’ 
formulas for certain polynomial sets. The last section is devoted 


to proving, under some more assumptions that the zeros of polynomials 


that are orthogonal with respect to convolution are real and simple. 


4.1 Convolution Orthogonality. Let w(x) be a non-decreasing 
function of bounded variation on (-~,~) such that all its moments 


fee) 
J tJay, y se051,..0, -exiistvand has infinitely*many "points ‘of 
—_—0O 


co 
increase; that w is a distribution function; and let {p, G3, 


be a polynomial set orthogonal with respect to jy, that is 


s 


(4.1) fata) By Cardv(e) = 1585 


and let L[f5a_ x] bevane (G,0)-- eranstorms of Q,() is the 
preimage under (L,a_) of P,(), then the orthogonality relation 
(4.1) is equivalent to 


(4.2) 2B SL As ie IGS 1 oy Fe 


Relation (4.2) expresses a kind of orthogonality for the polynomials 
Qo (*) 5Q, Gx) » ... . We shall refer to this type as "Orthogonality 


with respect to convolution" or simply "Convolution orthogonality". 
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To be specific, let (9, Go}, be a simple set of polynomials. We 
say that this polynomial set is orthogonal with respect to convolution 
if there exist a distribution function w(x) and a convolution * 
induced by an (L,a) transform that satisfy (4.2) for some sequence 
Aosryores of positive numbers. In other words a polynomial set 
{Q,G0}, is orthogonal with respect to convolution if and only if 
{LIQ,sa, x1}. is orthogonal in the classical sense. Note that 
one can drop the non-decreasing requirement of ww and still talk 
about an analogous concept. However, we shall always assume that 
y is non-decreasing. 

If we know a Rodrigues type formula for P00), then by 
virtue of the above correspondence and (3.6) we obtain a correspond- 
ing relation, for Q,™ involving BO » G@fscourses ,thiciwill be 


made clearer in the examples of the next section. 


4.2 Examples. As a first example we take {P00}, to be the 


Legendre polynomials of argument 1-2vx, that is 








= er st+-1T:]:% 
(4.3) es oF Cisittilive) 
Ce) (Pp) (1) 
and we take La3) ge he The ~ (i; at ) transform maps n to 


lieth yes t so that 


- (p) 
(4.4) SC GAG3 eats beta tGn Per wh Ewe 


p is not a negative integer. 


The Rice polynomials B, (5 3P 5? are gf o(-d,5+1,851,p.v), see [31], 


p. 287. These polynomials were introduced and studied by Rice [ 32) 
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imehicocase arelotioms4.2)¢ror .v. > di~ois 


(p)_ (@) 
m n re | 
(495) a eeu n,p,v)H, (-m,p,v) Dae eet 8, g/v(2i+1), We if! > 
with 
Liv 
We = taipe f sb aiees k ae : 
fe) 


The Rice polynomials H (EP .v) are not orthogonal in v unless 
Pe 1 or Wn. 4S. tollows reomese OfliLi. ginethe cases ££ = 1 or p, 
they reduce to special Jacobi polynomials. They are orthogonal in € 
fieand only iftiv.=815° +The«1f part’ folf@ows from,86\ef/ [1] and’ the 
only if part follows from Theorem 4.3, to be proved later. 

Thus the Rice polynomials for v >1,+-is an example of 
a polynomial set that is orthogonal with respect to convolution and 
not orthogonal, even in the generalized sense [33]. 

Rodrigues’ formula for the Legendre polynomials PO) 


Teme ot | ap. LoL) 
(4.6) P(x) = iG Oat) 


The BY operator corresponding to ans (p) ,/n: is -AnV-(p-1)A. 
Therefore, by (4.3), (4.4) and (4.6) we get 
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The Bateman's polynomials F C4) are H(z,1,1). In this case (4.8) 


reduces to the rather simple form 





, aes. | 72. p28 
(4.9) FC) = (3)1 (Azy) (6 j ) 


Our second example is the Meixner polynomials m,.(x; ,c) 
z 


ferined by C195] rvole -2,ap. 225). 


aan | 
(4.10) m (38,,¢) = (8), oF (-5s-%5B31-c ) 
They are related to the Laguerre polynomial Lago as 
(4.51) Lm, (n38+1, b+) 3(-1)"/n! ,x] = j! us” (ox) : 
since 
(8) (B41)*. 
(4.12) be) = pr Fy 58+) 3x0) 


The orthogonality relation (4.1) for the Laguerre polynomials is 


(eA oes ray 7% 


: j Q (bx) dx = b 
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As a final example we consider the transform (E2(-o a) 


that maps Bie) to ee . The Charlier polynomials oe are 
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defined by (see [34] 


: 4 j F ; . 
Meera aiatta 400) GeO) Ck a 
J fe kek 
saiG) 
The simple Laguerre polynomials Ei) are te GERCe. 2). 
Clearly 
(4.16) Le, (nya) 3(-1)"/(n!)” x] =ehonyal* eq: yad/2 1, (x/a) 
By (4.2), and (4.13) and (4.16) we obtain 
y gang te 
(4.17) Ne c,(n,a)c,(m,a)6, | = a 6, o/is : 
m,n=o 
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_omin ~ (-a)” ,mint2r 
we i. ) i ( mtr ) 
r=0 
Note that 
ee ee ee 2a 1/25 (5) 
i Ie nt (Bt1) ?* (-1)~ a Cpe) = ( /a)—ae 


leads to a relation of the type (4.17); with 67 ne of course, 
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depending on a; but is rather complicated. The 8 operator 
a 
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we get 


(4.19) e,(x,a) = Gy? ental “trea tc 62? 


RelLalioneaC4.0) 9°47) e(4 <0) 0iC45 9), C4ele) 2 (Aeby) cand 
(4.19) seem to be new. 

In contrast with the Rice polynomials, both the Charlier 
and Meixner polynomials are also orthogonal, in the ordinary sense 
(see [15], vol. 2, pp. 225-226). Furthermore, both the Charlier and 
Meixner polynomials are preimages of Laguerre polynomials under 
(L,a) transforms. This suggests the following problem. Given 
an orthogonal polynomial set {p,()}, » determine all the orthogonal 
sets £0, (x) 3, that are preimages of (p, (0) 35 under (L,a) 
transforms. The answer to this question in the case {P,() 3, are 
the Jacobi or Laguerre polynomials is given in Theorems. 4.1i.-and 4.2 


respectively. 


Theorem 4.1. The only orthogonal polynomial sets that are 
preimages of Jacobi polynomials ([31], p. 263) of argument 1-2vx 
under (L,a) transforms are 3Fy (a nty, x38) 58531). 

Proof. Let (QF be such a set. Thus Q 6) is 
of the form 
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(4.20) op OS ee oo ia ety bra Dh w emers 


where he ee tor 1k = Oy] acoae “Let 


Qty oe) = CAQxtB Cx) + CQ Ge) 


be the three term recurrence relation satisfied by QO, see [1 ] 
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and [32]. From (4.20) and the above three term recurrence relation 


we get 


—(n+1) (ntyt+k) (nt+¥+k-1) = (B_-kA_) (-ntk-1) (nty) (ntytk-1) 


2 2 -nt+k-1) /n - , 
C(-ntk) (-otk-1)/n + A (nty)A, j/a, 


It is obvious from the above equation that 


k=1 | 


2 
ermind ko tak tbicte =, a,.b,.c.. are constants 


> 


and where AO) is 4F3 (monty .x,158, 8, 8,51) and the rest follows 


by Theorem 4 of [1]. 

Theorem 4.2. ‘The only orthogonal polynomial sets that are 
preimages under (L,a) transforms of the Laguerre polynomials 
og OO sa. are the Meixner and Charlier polynomials. 

The proof is very similar to our proof of Theorem 4.1 and 
uses Theorem 3 of [1]. -See also [20]. 

Note that Theorems 4.1 and 4.2 are also extensions of 
Theorems 4 and 3 of [1] respectively. 

Now we come to preimages of the Hermite polynomials under 
(L,a) transforms. There is no orthogonal polynomial sets that 


belong to this class. In fact we have the more general result. 


Theorem 4.3. There is no orthogonal polynomial set of the 


type 


/ 
QQ & = } 


The proof follows immediately from the three term recurrence 
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relation. 

Theorem 4.3 tells us precisely that the preimages of a 
simple set of symmetric polynomials never constitutes an orthogonal, 
in the generalized sense, set. 

Now, it is clear that if we specify the orthogonal poly- 
nomial set tp (x) then the same question could be handled as we 
did for the Laguerre and Jacobi polynomials. However, the problem 
of determining all pairs of orthogonal polynomial sets 
(£2 (x) }o, {Q, (x) 5) that satisfy LIQ, tm) 34x] = ey for some 
(L,a) transform, seems to be much more general and remains to be 


investigated. 


4.3 Zeros of Polynomials Orthogonal with Respect to Convolution. In 
the present section we restrict ourselves to sequences f for which 
) ak remains bounded on (0,”). Thus definitions (331 and .(3.4) 
are equivalent. Furthermore we shall assume that $(x) is a Laplace 


transform of a non-negative function, that is 


(4 aah) o(x) = i} el e(t)dt fon exhe (0,0), with teat) > 0 
oO 


for tte (Ogee... 
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(4.22) Oe lesa) ete) 1 (t) de 
: Oo Oo 


An operator T is variation diminishing if V{Tf} < V{f} , 
where V{f} is the variation of f defined as the number of sign 


changes of the function as x varies across its domain [18]. 
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k 


Lemna’. "hét £(x) = ) a < be uniformly convergent on 
Oo 
fesar, thenfon?? [0 a], °vtr] = Vita, }], where Vita], is the number 


of sign changes in the sequence AysAyseee 5 


Lemma 2. If f(x) = f et o(t)dt, then on [0,~), 


0 
Vitle<.V[el- 


Lemma 1 is due to Cheney and Sharma [12]. It follows easily 
from Descartes rule of signs. Lemma 2 follows from Theorems 7.1 and 
9.la of [18] on pages 97 and 103, respectively. 


We now come to the main result in this section. 


Theorem 4.4. Suppose that iy Gr 2 is a polynomial 
set orthogonal on a subset of [0,~”) and Q, Cx) = L tp, (#) 34,0] 


with (x) = ) ox satisfying (4.21). Then Ci) has real 


and simple zeros. 


Proof. Lemmas 1 and 2, and (4.22) imply that the number 
of sign changes of Dale) in,  [0.@) | sis=at most Vi{Q,(n) }_ I. 
On the other hand Be) has j changes of sign in [0,”). Thus, 
j sign changes of On 2) muctrocctih at, x =. Oo... tand the 


result follows since a) is a polynomial of degree jj. 


Corollary 1.) the zeros of the Rice polynomials are real 


and simple for p > 0. 
Proof. One can argue that 
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BO ee ure: x] a Gee 3 


P(x) being the Legendre polynomial of order j. For v>0O, _ the 
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Legendre polynomials of argument 1-2vx are orthogonal on [0, a 
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The corresponding $(x) is (1+x) ? which is the Laplace transform 
-£ a 
of e tP eens Therefore HLCP ey) has real and simple zeros 
fone wae UAr andr: p> -0. Let 
j ; j 
- 2 
Ce ( J) I (x-6 


: ans Vino ete 
aa anaes! k,Psv 


ce Ce.D Vv) = 


where the 6's are distinct. Clearly the 6's are real polynomials 


in yaa Therefore Seis ie, also has real and simple zeros. 


Corollary 2. If {Q,G0}, are orthogonal with respect 
to convolution and the points of increase of the distribution w of 
(4.2) are contained in (0,%), then the zeros of Cee) are real 


and simple. 
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CHAPTER V 


MISCELLANEOUS RESULTS 


This chapter contains several miscellaneous results 
related to the subject of this thesis. 

Recall that the theory developed in Chapter LIL enables 
us to handle difference equations of the type 


ml 


Ch) ) aT “£(n) = g(n), ay24 


Dee) are constants, 
om m 
k=o 


er’ 
for operators [ of a particular type. Consequently one would 

like to have a way of recognizing difference equations of the 

above type. In other words one needs some operational formulas 

to expand powers of [| in terms of the shift operator E and 

its inverse ay in order to be able to write a given difference 
equation in the form (5.1), if possible. So, we devote the first 
section, 5.1, to expansion formula of (AxV)” in terms of powers 

Of. EE. and mee and illustrate how to obtain such formulas for more 
general ['s. Osipov ([27], [28]) studied the problem of expanding 
the Bessel operator DxD and its generalization shea in powers 
of D. Carlitz [11] independently gave a similar expansion for the 
Bessel operator. 


In §5.2 we introduce a q-analogue of the Bessel operator 


DxD, namely the operator a Zar where 
(5.2) ie) = {£(qx)-f(x)J/x(g-1) 
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As we shall see these considerations will lead to a q-analogue of the 
Laplace-Carson transform. Some applications related to Heine series 
and to q-difference equations will be mentioned. It is clear that the 
differentiation operator D is the limiting case of Pa SS08d ua, 
Finally in §5.3 we shall indicate some possible extensions 


to multi-dimensional cases. In particilar we shall indicate a method 


of solving some mixed equations. 


5.1 Expansion Formulas. Recall that 8 = AxV satisfies 


BE(x) = Cxtl)£(xtl) — (2etl) ix) - xi(e-1+)”’. 


It is clear that 


(5.3) eg) CD ae Gye, 
k=-n : 
that is 
n = k+n 
(5.4) ®2£(x)m= is) .(-1) ane Ope) e, 
==—N 9 


where the a pO) 'S are polynomials in x of degree at most n. 
Nn; 


The relation 


Coe) 


oO t-Ih 


fea ERA AN 
(1) ()() = ae > 


for non-negative integers r and s, r>s, follows easily from 
the binomial theorem. 
Let f(m) = > in (5.4), for positive integers mn, 


to get 
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ap tise kn mtk 
(5.6) C= eee L GLY a, eT). « 


d ; sure Banal j 
Multiplying (5.6) by ‘Bye < J ) and summing over all j, by 


(5.5), we obtain 


Cee lear Gilaean! lame) Cee) Clu: 


eG Saati ya here aig ect TD 


In the left hand side of (5.7), j runs from min {n,m+2} to mtn. 


Replace j by jtmt2 in (5.7) in order to have 


mae sy gtn+2 ,jtory j+m 
(528) ay gm n: ) (-1)~ ¢ cages »¢ pe 
j 
Now, for fixed n and 2, the polynomials aay a ge /ne and 
) a ORES ae Oe. agree at all positive integers, hence 
j 


they are identical, that is 


eo morn Gesell) 2) C1) 1 


where 


C) Pex (Rae Agr) fy 


Clearly the above procedure can be duplicated for any 


operator [ satisfying 


m = = 
rq) = eu D> Pe ek a 


Indeed let 
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(5.10) ete GL ee EE (x) 
tae nk 
mM 
For f£(m) = (> we get 
1 ae - k+n m+k 
(5.11) ese, yee, yy Gi) = eA: ak.) 


k=—-@ 


Multiply (5.11) by Guy 4d J 


and sum over j to get, by (5.5) 


Node aney m 
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Weysee thatel {a caeseks, 1.6.$eehi= 8 , formulary (5.12) 
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reduces to (5.8). In the case [. = of Chapter III we have 
c. = Cytk) a One 


and 


a ae rae) 


(5213) ang =n! ) (-1 


en me 


Clearly (5.13) reduces to (5.8) when y = 0. 


Now we go back to (5.9) to write a, yg) in a more 
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compact form. Clearly for k > 0 we have 
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and (5.14) imply 


nt+k ae 
Baio ye 9)eona, (=)18j ane coe Gahliv wal WRATH) Cay be1Qalone ses 
2 j =O ay 
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ils 
€5.16) aK) * tae An { (xe &™ 32 


Kee Oe. nk 
Furthermore, we have 


n-k 


(5417) a, () = y 
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Relation (5.17) may be written as 
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Kea OSL ss, 2 


Going back to (5.3) and using aor B gn we obtain the 


recurrence relation 
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2 
is A. b = (xt1) +(2> +h x--1) 
( ) Fea) . : ARONA hy oth el 2 Bere 


From (5.9) and (5.16) one can easily see that bt cor) also satisfies 
Oo. 245) lL hus, Lb (x) and ib (x-k) must be identical since 
n,-k n,k 


bY a = an = 1. Consequently 


) 


(x+1) 
(5.22) a, go) = Taye a YP 1 Rha ite Dee ee 
and 
(5.23) Seer Bethe Ce aed (ees he et Oh lS in 


This leads to the curious identity 


RL SORE ag sees 39) a bode 


5 


Cye24) ieee Wen Ta wea ee = { (#1) 37 (nk) 1 


5.2 q-Analogue 


5.2.1 q-Analogue of the Bessel Operator. The q-difference operator 


De is defined by 
Df) = {£(x)-f(qx) }/x(1-q) 


It is clear that, for differentiable functions, the differentiation 
operator — is the limiting case of Des as q-2>kl1. We define a 


q-analogue of the Bessel operator by means 


£(x)-2£(qx)+£(q-x) 


(5.25) Buf(x) = (D-xD )f(x) = 
: ied Cea 


We begin by introducing some notations. q-factorials 


[n]! and q-binomial coefficients [a are defined by 





64 


(5.26) [n]! = (qsn)/(i-q)” ’ , 

where 
n 0 

(5.27) (qjn) = H (1-q) , n>O and (q30)=1, 
i=l 

and 
is in bees Ee n 

(5.28) Mee teyacei i teed 


Moreover, by [a] we shall always mean (cp) iqel)e We shall 
assume O <q <1. 

There are two q-anailogues of the exponential function. 
Indeed the exponential function e* is the limit as q>1 of both 


FOSS, and oe oe where 





e Cie gel pe 
(5.29) Soe meeie (lox) gudeel) x /{qsi) 
4 =O N=0 
and 
co . 2) e ial 2: ne: ; 
(5.30) E (x)= 0 (iteql) = J GSD! siy(qsiy 
q 1=0 j=0 
We shall also use [17] 
Dare -v 
= Pad Ce Glee iy 
ae a aa 
caret qt n=o A 


In particular if v=+N,~-N is a positive integer, formula (5.31) 


becomes 
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| ae N . 
eer Gt) Gia ata = yea 
nh=o ¢ 
and 
N co -N 
(5.33) 4+ = ) alae ye 


(aid eee Cathay n=o © ~ 


respectively. 
We are interested n finding a transform, oe say, and 


a convolution produce * such that 
(5.34) rg ites! = Ei lfsxll lesa] , 
holas.and Jf jmaps B. to D , ‘thateis 
q q q 
HAs Li Bhs: xhistpebel( tf: 
( ) a - ] qi x] 


In this setting, one might expect that 


(5736) L (eh /Eo] 3x] = ds 
x" sat n n-L 
This is so since BG Tape = [n] (enn S while tips = [n]x . The 


above interlation between *, op tT and DF resembles the 
interlation between the convolution of Chapter III, the (L,a) 
transform, the operator BO and the differentiation operator D 
respectively. 


We recall that the q-analogue of the Laplace transform 


(see [17]) 


2 n 


= al - ares q a 2 
dts Bsa ae #(—a) a CR Ae 


. 4 
— * te 


ae 


. GS ial = ae > 7 : 
eee ammcarreles SelGes 21) Mo sugotkaesp ag? ss aes Cina 
= : f a } 


_— lege * 





satisfies 


n —T) = 
Gees aa) Crem eth 2), 


which suggests that we consider, as the q-analogue of the L 


transform in Chapter II, the transform =f defined as 


(5.37) L [fx] =E(-q) J @ £29 1455) 
q q fs 


It is easy to see that (5.34) as well as 
52 Pe Dee x= tol esl ort 30 ; 
€5238) at P xiest a a | att ae 


aAre@walid under (5.37) 


Furthermore we have 


Theorem 5.2.5. The relations 


n tps n : 
(5.39) va Pts] = eee iy ed 
and 

dpa a di n y 
(5.40) Gita £(t)igx P= ty) a! 


apd =(Dea5) holds. 


Proof. We have 


nt1 


Ea(-d) a 
q xq 
624 pare bet 2 Aza ) 3 


Bk Bye es, ete Rie ae 
xD {xl lfsx]} iq 





li o1 8 


q 
path (q3n) 


ee eG } te 


Th 
oar: 


co 


BE (-4) L TER OEE yee = Li[et(e)sx] 


I 
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and (5.39) follows for n= 1, and by induction for all n. 
Formula (5.40) is obvious and (5.35) follows from (5.38) and (5.39). 


Let LD be the inverse of me We have (see [17]) 


SB Oe add Ce) 2 ; 
Peace = ahd) ) 2 fsa) 1-4) 
Oo 3 


Cac41) 
Note that (5.37) as well as (5.41) are valid at least for 
functions that are regular at the origin. 
We are now in a position to introduce a convolution 


product * so that 
(5.42) [ ltxesx] = pelt stle ber) 
Or, in other words 
(f*g) (x) = fail bebe ee aes, 
q q q 
By (5.37) and (5.41) we therefore must have 


4 (-1)4 quer: VP ce 
(5.43) (f*g) (x) Glee (=a) ) Sie 85 ere he 
q 4 (q34) sty ie 


2 ae £ (xq) g(xqhts) 
(q;m) (q3n) 


Now having the definition (5.43) at hand we can verify (5.42) 
directly.turhe set of all functions that are regular at the origin 
forms a commutative ring, under addition and convolution (5.43). 
This ring has no zero divisors. One can easily construct an 
operational calculus based on the convolution (5.43), as we did 


in Chapter III. The technique is very similar to that of Chapter III 
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and we shall leave it out. 


9.2.2 Application to Functional Equations. First of all, the linear 
q-difference equation with constant coefficients is transformed 

under te to an algebraic equation and its treatment is identical 
with our treatment of the linear difference equation with constant 
coefficients in Chapter II. 


Our second application is to q-difference equations of the 


type 
n : 
(5.44) ) a. (xlax) FCs) = C(x), 
ae J q 
j=o 
= | 
where:the a's are constants. Let f(x) = on [F;x] and 


g(x) = Eales! Equation (5.44) is therefore equivalent to, by 


(5.34) 


nh - 
(Gia te fCO. = 2(t) 23 
oO 


and hence 


The limiting case, as q7*1, of (5.44) can be treated similarly 
using the Carson-Laplace transform. 

Our iad application is to q-diffierence — differential 
equations with constant coefficients in er and x a Such an 


equation, say 
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y-£(x) + )} bd £(x) = g(x) 
a 2 jar 4 


is transformed under FT to the much simpler equation 


f m ; 
A a, (x <)* ECR et (} y. x J) F(x) 
i=30 jal 


M15 


tr 


nj 
vibes 
a 


pet eo Lie Gy cage, 
j= = 


EO 


F(x) and G(x) being the be transform of f£ and g_ respectively. 
As an example of the above type consider the functional 


equation 


2 


(Goya Gee eldx) Cs) = x Ag-1) ) een (ta) ' 


whose Fs transform is 


Se oa + x/2 
x a 


xF" (x) + 


where F(x) is as before. Thus 


FC) =) E(O) + 


Nie 
™ 
w 


or 


Exe 60) eae 


5.2.3 Applications to Basic Hypergeometric Functions. A basic 


hypergeometric function a is defined by the Heine series 
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2 (1-a,) +++ Cea ) 


a mpg = mi Mera 7. 
mon 621 2°71 24,3P] 92 9. 3x) a (1=b ) ...(1-b ) x GEC) 
r=o hae Dot 


These functions are q-analogues of the generalized hypergeometric 


functions ia 
mn 


Meller [22] and [23] used the following method to obtain 


identities among special functions.) Let qT) and T, be two given 


operators and suppose that one can find an invertible operator U 


satisfying UT, = TAU and Ul = 13 where 1 is the constant function 

f(x) = 1. By induction we have ur = Tus Toes late alneretore 
=n 6, 4 —n ,. -n 

(5.45) vd. ae) De Can De) Ul ate) at Ti 5 


In general (5.45) and 


(5.46) eat, D1 0 Ge wre - 


express identities. If we take for Tis Ty; q-difference operators 


we get identities involving basic hypergeometric series. In par- 


n 


ticular, put a = aay that is 
cs 1 eee 
) rer Toa 
fe) 
=18) Sal . 
< = = ; th 
Case: f. Take qT) on and T, yaa Dt In is 


n D 
case it is easy to see that Ux? = (1-q) ((1=q) /Ci-q) ae es 


so that U may be defined by 
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' esky itneb? sesege 
itvinwal eataliwebhk 30g ae 
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fee 


ee Bj _4(j+1)/2 
(5.47) f(x) = eG) F (qt ge 


a F(xqI(1-q)) 
a2 (433) 





and one can check directly to see that the operator U defined 


by (5.47) satisfies all the requirements. Furthermore, we have 


a On A Biss 5 
U wie = ELCd” B) ; A £(=4 


eae ds I=q 
a Ty 2 
T-a + = a) and ‘T-a 1= 9? 16-3 8t13 Cl-q) ax) 
From (5.45) and (5.46) we get 
= Saltese 
B+1 Bj 4q a 
e€ is ; = = etis 
ACIS ho (-q°) Casiy 8q S24) = 996-3 8t1s%) 
and 
se eri ' 
Bt1 (Uae) j 
= E (- An ~;B+13x 
en) = Ecc ”) h Caray 081 6 38413x0") 
é ee _ 6-1 B+1 
Case II. Take tN Dy and T., x ee ; 


# 0,-1,-2,... . As in Case I we first establish 


(1-4) 9 (asm) 
= a OK 
(l-d) py 


Th n 


Ux 


and try to guess the definition of U. Note that the operator U 
of Case I maps x to ((1- Cle )_)/CC1-qa) ) x" arene ge 

a map q eee q ren a q 
maps x to (quia tied): Thus it is conceivable that the 
composition of these operators is the required operator. This is 


indeed so and we set 


if 7 ? , iar 
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= y | : : i be s . 
{ 2 { } > es 
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me ahd Lee B+j-1 q 
UE(x) aw FEE ecg), 
and hence 
, B+1 
E (- oo = 
rege so eas ) aaiiad cr £ (xq) 
(1-q) 44 x (q34) 
In this case, 
uy oy) 
T-a 1 = e,((1-a)ax) and T,-a US 191643 Bt1; (1-q) ax) , 
and (5.45) and (5.46) imply 
B+1 
ES (-@) cata = gy 8 "Ca 
eee ee ee een) 
(1-4) 94 O : 
and 
Bars 
Bed gee er (gy om, (Clq) ; 
q q ie B+jtl reer See 2 


It is clear that we can repeat this process and get many 
more identities. For more results of this type regarding ordinary 


hypergeometric functions see [22] and [23]. 


5.2.4 Another q-analogue. The 7 transform of the previous section, 
being a function to function transform, might not be a proper analogue 
of any (L,a) transform. However, one can define such a q-analogue 
by 


is n 


(5.48) Al [f3a,x] = ) (=x) -a(a-1) 
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co 
where a stands for the sequence ta to and (x) = , 
i 


Oo m18 


n —n n —n 

a $(xq ) we mean, of course, DF @ evaluated at xq . In the pre~ 
sent section we shall restrict ourselves to q>1, and q=1 will 
be a limiting case. If $(x) has a positive radius of convergence 


then 


Ss 
a x® ; higibe es 
fe) n=o 


n(1-s) 


Cyc 47) A {fsa,x] = £(n)q ’ 


i 1 8 


Ss 


for x belonging to some neighborhood of the origin. We shall assume 
that o tee nes O. be. ae, ew nsorder to make ih ont-to-one. Further- 


more we suppose a. = 1. Let 


j(j-1)/2,n) 


5 50 as = od : = (-I Jo (n)/a, : 
(5.50) yo q ; Bm) i) j a, 
and 

(5.51) ca 0 y Sand Coad =~ M4! o pee 10 


It is easy to see that 


n 
(5.52) iin amen mee CE) 10 Car een (e051 acs) us, 
quel 3 j 
where 
(5.53) b.(f£) = (quea)® £(2) 
J L=0 
Formula (5.53) may also be written as 
Gis) by) emg aac 


eee) 


and symbolicially, by (5.32), as 
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b,(E) = qIG-Dl2 (1) (H-q)...(E-q! 4) £(0) 


Now substituting for f(n) from (5.52) in (5.49) and using 
(5.32) we get 


oO 


Ss 
dee) uci Deh) 


i] 


(5755) A [sa,x] 


= rupees day ASD)! (aan cay. yea eco) , 


which may be compared with (3.4). We shall adopt (5.55) as our 
definition of the WPS cogs or oe tranfornm. 


We point out here that (5.55) may also be written formally as 


foo] co 


: = pee ot (i) 7 2 mp ety ie 
eos oo be 9) q i (ty) ne [ r lq Caer 
co 1s 
={) (7) ——BO_____jyz , 
L=0 (1-xE) (q-xE)...(q'-xE) ° 


by (5.33). Therefore we get the formal relation 
A (f:o.x)] = A [£;(-1)" -xE] a 
q > 3 q 3 3 fo) 3 


which in turn may be compared with a similar formula in Chapter III. 


fab Wy . : 
The case pee (-1) is of special interest. We have, 


cox | x| —e1, 


2H n 
(5.56) A [£3(-1)",x] = ) (qx) EGO 
: n=o (1+x) (qtx)...(q +x) 
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n foe) a 
Ate Soe, xb =e (abs) Y Sry 
: r=o [r]! ((l-q “)xtl)_ 
beet) a. = 130 n(n-1) /2 
r ale ae » we get 
(n-1) /2 ©o 
AU Sc ae ete) 
es [n]! cy x] Bota aD iar Tae 
On the other hand, corresponding to (3.36) we have oe Caleta 
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A, te eee Pes ) eo) ; -r(r-1)/2 £(x) 
r=0 yg 
co ~r(r-1)/2 
= J xX 2 (87) ey) 
a (itxq eg ae 


Let us define the operators A and B . that are 
O54 5G 


q-analogoues of AY and BY Om Chapter Lit.eby 
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are obvious. 


Theorem 5.2. We have 


(5.61) Bsr Reais = x) a erate 
(5.62) NL Meth] o 2 eee) 
q aq x 
and 
(5.63) AULA, gtstee = Le iG ae 


where F(x) = SES eae 


Proof. Relation (5:60) clearly follows from (5.55) and (5.52). 
The definitions (5.57) and (5.58) tell us precisely that 
ong. = Cra ae and ES aca = Been Cs ey, Peat? respec- 
tively. This obviously proves (5.62) and (5.63). 


Thus a linear equation with constant coefficients in A 
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is transformed under ee to an algebraic equation. 
We define the convolution product by 
m 
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eee Alma) Veh [nal (atl i+ (hen ein 1 Ven) elen] f(s) ts 


In the rest of this section we shall denote A fes(-D" x] 
by A [fsx] or F(x). Recall that A ffs(-D" x] is defined by 
(o250)0 > Glearly. sas "q 2.14 A lésxl, q A and aq? A{-n]¥ reduce 


Pou L{ taxi. A and 6 of Chapter @l,, Ltsis easy to deduce 


(5.66) A [é(ntl) 3x] = F(qx) + {F(qx)-FQ)}/qx , 
and 
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(5.67) esc e 9 5 ae i ies 


Proms 3.50) . 


Theorem 5.3. The following formulas 
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and (5.68) follows. Clearly we have 


a Mu (xq)” 
(1+qx) ) ers ete ee £(n) 
of (xt) Ga (aeq ) 


i 


F(qx) 


n n 
Cash , xq) __ Geta) £(n) 
O- Gxt lk< (ckg «) 


iH 


Citas yee {xP Gx) +A La" £(n) 5x] } > 


which proves (5.69). Formula (5.70), for k=1, follows from (5.69) 
and can be established for general k by easy induction. 


k 
We now proceed to evaluate the Ea transforms of E and 


-k 
Pe owe. being the shaft operator, By Ltenating 165 26/)esve, pet 
-k 1 k 
° A ¥ ] CS STE ea . 2 y 
(5.71) ae SS Pam eye 


k 


To evaluate Peele let hts a eit Datei. fa Bier {roar t ee 


f OO... 2 anus we have 


ie 5 de a 
k-1 j k 
(xq)" X -k 
F(x) = ) £3) =e — + 2 Gag) 
- (xt1) 4, (etl), 
or 

k=-1 j 
~ -k k. -k J 


-k 
(x+q ) 541 


The. (A RL} transform transforms a linear difference equa- 
tion with constant coefficients in q A to an algebraic equation. Note 
that any such equation is also a linear q-difference equation with 
constant coefficients. This is so, since by the change of variables 

-n 


q = x, and £(n) @%e(x)- q “AE becomes (q-1) a etx). 4 1ne 


inverse transform transforms a linear equation with constant coefficients 
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(5.73) ; ay {xD esa PCxyesrG(x)- /, 
k=o q 


j 
to {) ain} “}£(n) =Je(n), with £(n) = A [F:n] . and 


g(n) = nv Gen]. Consequently the solution to (5.73) is 


Av {e3n] 
F(x) = A | —4 x | oe 
qiset 
i: 
) a, [a] 
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Our final result in this section is to establish expansion 
formulas of powers of a general AL operator in terms of powers of 
> 
the shift operator E. This is a q-analogue of our investigations 


in 85.1. On the other hand the results of §5.1 may be looked at as 


bphetlimiting.case.,q 271 of our present result., Let 


k 
(5.74) Ay (des) CI 2, 8m) EE(n) 
of k Qk 4 
in particadtars Lf 1 (n= ee we get the relation 


Cae: tess Ce Or i) = a (n) @,(ntk) , 
sheet BE oa L & yk i‘ 


or its equivalent 
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Summing the above equations for different j, with mtr 7) Sart 
and using (5.32) we get 
ipeooas ghia pad oenel) 12 n 


[oe tad 


Car 75) om) = : (-1) : 
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ntr 


O(h-4-1) fe 


x Ce Ce eee Cj-g 4 


Clearly (5.10) and (5:12) are the limiting case q->1 of (5.74) 


ence US e753) 


5.3 Multidimensional (L,a) Transforms. A two dimensional (L,a) 


transform, say L[f;a »x,y] can be defined as 


m,n 
,mtn 
(5.76) EEC eee y ee ee Lao LES) Re aa 
m,n Meet. igen m,n 
99 ox oy 
on double sequences {f Tee » where (x,y) is the power series 


m,n (0,0) 


(5.77) o(x,y) = ) Sar ae ee 


We shall use the notation 


i 


A, £(a,n) f(mtl,n) - f£(m,n) 


£(m,n+l) = £(m,n) 5 


A,f£(m,n) 


and similar, notation for, Eb. , V andeV « 6 
uf 2 us 2 


The convolution product is defined by 
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uf ay +b+c+ 
DoE ty gions aad et Na a ead aan 
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C7) 8) O57) AF) £(a,b) Ce .d) 
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x A, A 9 £00) 2A a 2¢(0,0) 
Note that (5.76) may be written, formally, as 
(5.81) L[fsa,x,y]= ] (-x)"(-y)" a, Ay A> £(0,0) , 
0,0 
where a now stands for the double sequence {a We pa 


Mestl (0,0) 


We shall restrict ourselves to the case $(x) = Ctra 


Ena cis as = (-1) 


on the same lines. 
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a ap) and the general case can be developed 
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3 Pe) 
5.84 = —_ > — 
( ) Q= 1+ x se + y ay 


One can show that 


(5.85) L (mbnt1) (myn) 3(-1)™" (™") x,y] = (itxty) F(x), 
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HT} 
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In this case we have two finite difference analogues of the 
Bessel operator. Indeed A, (av tnV,) £(m,n) and A, (mv, t+nVv,) £(m,n) 


are such operators since 
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Furthermore we can prove the following formulas 
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k : . k 
oie A I (- V = i 
(5.93) ia (-rt+m 1) i (mV tnVi +r) A, ¢ 
r=0 r=o 
fe k 
and two more formulas similar to (5.92) and (5.93) where ‘Ay and <Ge 
k Ox 
k f) : 
are replaced by A, and ai aS In particular one can conclude 
: oy 
3° 3° 
that the Laplace operator sip Say is the transform under 
ox oy 
n /|mtn js 


‘CA aT He stoe (hat 5) (av, +nV,-1) (mv +nV,). 


The general eas transforms may provide a simple method 
to solve difference equations in two variables. However, the present 
topic is far from being complete and will have to be investigated, 
in more details, in future. 

We would like to point out that a combination of (L,a) 
and Carson-Laplace transforms may provide a very effective technique 
to solve functional equations. For example the integro-differential 


difference equation 


df(m,n) rs d 


n 
(5.94) £(mt1,n)-£(m,n) + an ae f k(m,n-t)£(m,t)dt = g(m,n) , 
Oo 


is transformed under T, 


fore} co n 
-1 i Rat x Bie td oe 4. 


(5-95) beer yay. 
: = Gia 


to 


B(x y) -F(05y) oe F(x,y)-F(x,0) 
x y, 


+ K(x) EF Oey) = GG yee 


where F, G and K are the T-transforms of f, g and k_ respec- 


tively. Therefore 
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_ xyG(x,y)+xF(x,0)+tyF(0,y) 


5.96 F ’ - 
(5.96) (x,y) xyK(x,y) tytx 


Now in order to get f(m,n) form (5.96), one has to be able to invert 
the transform (5.95). Such inversion does not present any problem 
since the inverse is a composition of the inverses of the Carson- 
Laplace transform and the L_ transform. 

In general one might be able to handle functional equations 
in several variables by using a different transform for each variable. 
The extension of (L,a) transforms from two to higher dimensions is 


rather obvious from (5.56). 
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